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Abstract: This study provides a general framework to analyze the effects on correlation 
radiometers of a generic quantization scheme and sampling process. It reviews, unifies and 
expands several previous works that focused on these effects separately. In addition, it 
provides a general theoretical background that allows analyzing any digitization scheme 
including any number of quantization levels, irregular quantization steps, gain 
compression, clipping, jitter and skew effects of the sampling period. 
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1. Introduction 

Microwave radiometry is today a mature technology that was first used in radio-astronomy in 
the 1930s [1]. Since then, a large number of microwave radiometers have been developed for remote 
sensing applications to measure a wide range of natural phenomena (for examples, see [2-5]). 
Continuous technological evolution has given these systems new capabilities and features. One of the 
most relevant new technologies in the 1960s [6,7] was the digitization of the signals and the 
capabilities that emerge from specialized processing platforms such as today powerful an omnipresent 
Digital Signal Processors (DSP) and Field Programmable Gate Arrays (FPGA). Although digitization 
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provides versatility, re-configurability and other advantages for microwave radiometry, it also has side 
effects that must be carefully analyzed. 

Digitization effects can be separated into quantization and sampling effects. Mainly, the effects 
related with the quantization of the input signal (thermal noise) imply the loss of its statistical 
properties due to the non-linear quantization process. Consequently, it is not possible to apply the 
well-known Gaussian statistical relationships to the quantified signal. This effect has a large impact 
when limited quantization levels are considered, and it can be mitigated by increasing the number of 
quantization levels. Non-linear effects studies on Gaussian signals started with the early analysis of the 
spectrum of clipped noise by Van Vleck and Middleton [8]. In the late 1950s, Price published a work 
focusing on the relationship between the ideal correlation between two random signals with Gaussian 
probability density function (pdf), and the correlation measured after a non-linear manipulation of 
these random signals [9]. This relationship is now used to study the effects of arbitrary quantization 
schemes on the correlation of two signals. 

Sampling has an impact on the correlation of the two input signals by creating spectral replicas. 
Additional noise can be added to the mean correlation value due to spectra replication depending on 
the ratio between the sampling frequency and the signal's bandwidth. Moreover, the sampling period 
and its inaccuracies (skew and jitter in the sampling periods) of the Analog to Digital Converter (ADC) 
can affect and distort the sampled signal and so the cross-correlation value. 

To compare different sampling and quantization schemes an equivalent integration time is defined 
as the one required to obtain the same resolution as in the ideal (analog) correlation. Hagen and 
Farley [10] conducted important work on this topic, where the effective time was defined and some 
easy-to-calculate digitization configurations were analyzed in depth. 

More recent works have extended the initial digital radiometer proposed by Weinreb in 1961 [6] for 
auto-correlation spectrometers in the radio-astronomy field to include total power radiometers, 
polarimeters [11], and digitization impact on interferometric radiometers [12-15]. 

This work provides a general context to analyze the digitization effects on the cross-correlation of 
two Gaussian random signals in the most general case including any number of quantization levels, 
irregular quantization steps, gain compression, clipping, bandwidth, sampling frequency, and the skew 
and jitter inaccuracies of the sampling periods. 

In Section 2 a general analysis of the correlation over non-linear functions applied to Gaussian 
noise is given. In Section 3 quantization effects are analyzed as a particular case of the results of the 
previous section. Section 4 analyses the sampling process. Section 5 evaluates the correlation variance 
due to the digitization. Finally, Section 6 presents the main conclusions of this study. 

2. Non- Linearity Impact on the Correlation of Two Gaussian Random Signals 

The input signal of a microwave radiometer comes from natural thermal emission. Let us consider 
here two independent Gaussian random processes with a certain cross-correlation coefficient 
characterized by their probability density function (pdf). Furthermore, both input signals are assumed 
to be statistically identical and stationary, so the properties of the signals at a given time (To) are 
independent of To. Furthermore, both signals fulfill that their statistical properties can be deduced from 
a single, sufficiently long sample of the process (realization), i.e., they are ergodic processes [16]. 
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Thus, the variance of each series is constant and the covariance between elements depends only on 
their time separation. Finally, let x and y be input signals of a radiometer fulfilling the previous 
requirements. Usually, these random variables have zero mean (u^ = 0), except for offsets, and are 
considered to have a standard deviation different from one . In this case, the joint probability of density 
is defined as in Equation (1): 



f{x,y)= e 2 Wy^-dy) (J) 

2UO X OyJl - Ply 



where: 



• a x y are the standard deviation of x and y signals, and 

• p xy is the Pearson correlation coefficient and it is defined as the ratio between the covariance 
of x and y and the geometric mean of the product of their variances (Equation (2)): 

(xy*) 

PXV J(xx*)(yy*) (2) 

where the brackets (•) indicate a statistical average, over a time sequence with identical statistics 
(ergodic properties). Price's theorem [9] relates the correlation coefficient of x and y after (R xy ) and 
before (p xy ) the non-linear functions (S^O) and S^O)) to the random variables (Equation (3)): 

— ■■- a*a*(gZ(x)g$(y)) (3) 



d Pxy 



where: 



• R xy is the correlation coefficient of x and y after the non-linear functions #i(0 and fl^G) are 
applied to x and y, respectively 

• p xy is the ideal correlation coefficient of x and y (before the non-linear functions S'iCO 
and 0 2 Q) 

• o xy are the standard deviation of x and y signals, (for the sake of clarity, from now 
on a xy = 1), and 

• k denotes the k th derivative of each function. 

From Equation (3), it is possible to retrieve p xy , the ideal correlation between x and y, by 
^-derivation and further correlating the g^x) and g 2 (y) signals and then ^-integrating with respect to 
the p xy variable from -1 to +1. These operations can be defined as a relationship between the ideal and 
the non-linear correlations. It is defined here as the function q [•]: 

Pxy = ^[Rxy] (4) 

Considering the correlation of two signals at a delay different form zero (x — > jc(t = 0) and 
y — > y(t = t)) the effect of the hardware response can be included as depending on t, then 
Equation (2) becomes: 

p xy (r) = (xy*) = p X y(0)r xy (T) (5) 



Sensors 2011, 11 



6069 



where: 

• Pry(0) is the Pearson correlation coefficient (Equation (2)) at t = 0, and 

• r^ y (t) = / B ^ By / 0 Hx(f)Hy(f) ei ^df i s the so-called Fringe- Washing Function 
(FWF), as defined in [17], where B A and B v are the equivalent noise bandwidths, fo is the 
central frequency, and H x (f) and H y (/) are the normalized frequency responses of the x 
and y channels, respectively. Note that, if the channels' frequency response are equal, with 
rectangular shape of bandwidth B, and centered frequency fo, the fringe-washing function 
reduces to the well-known formula ? X y(. T ) = smc(B xy r). 

Equation (5) depends on the correlation coefficient and the hardware features. Usually, to avoid an 
underestimation of the correlation coefficient, Equation (4) has to be compensated by the existing 
delay between x and y as it follows: 



Pxy(0) 



l xy 



CO 



(6) 



Assuming that the x and y are stationary random processes (thermal noise radiation) that fulfill the 
ergodic property, then Equations (2) and (5) can be calculated in practice using the cross -correlation 
technique (multiplication and time averaging). 

3. Quantization Impact 

This section particularizes the previous analysis to a generic ADC function. The degree of 
non-linearity of an ADC function depends on the number of quantization levels and the ADC span 
window (Vadc =X m — X 0 ), going from 2 levels (1 bit) the most non-linear scheme, up to an infinite 
number of levels with an infinite Vadc, which can be considered linear. As explained before, 
quantization does not take into account the sampling, so that this analysis considers an infinite 
frequency sampling (F s — > co). Sampling effects will be added in Section 5. 

Figure 1. Generic analog-to-digital transfer function. The function plotted has compression 
gain and different steps in analog and digitized domains. 
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Equation (7) shows a general quantization function (Figure 1): 



Af-l 



g t (x) =^A i m u(x-X m ) (7) 



m=0 



where: 



• x is the value of the signal in the analog domain (continuous in time and values), 

• X m and X m+1 are two consecutive threshold values of the input signal in the X space (see 
Figure 1), where X is the set of all the input threshold values in the digital domain 
(X = {X m }, m = 0... M — 1). Note that, in the most general case the distance between two 
consecutive steps does not have to be constant (X m — X rn _ 1 X m+1 — X m ), 

• u(x — X m ), is the Heaviside function (step function) centered at X m , 

• A l m is the quantization step value between the (X^^X^J interval. Note that, in the most 
general case the distance between two consecutive steps does not have to be constant either 
(Ajn * A m+1 ), and 

• the values X 0 and X M _ X define the lower and upper bounds of the quantization window. 
Again, in a general case, the lower bound can be different from the upper bound 

(Xm ^ — X M _ X ). 

The first derivative function of Equation (7) is given by: 

M-l 

dgtw 



= ^A^5(x-X m ) (8) 



dx 

m=0 

where 8 is the Kronecker's delta: 8 = 1 if x — X m and 8 = 0 if x X m . 

Substituting the results obtained in Equation (8) in Equation (3), with different g (•) functions for 
x and y, and considering the first derivative function (k = 1), the relationship between the correlation of 
two signals before and after the quantization process is obtained in Equation (9): 



dR xy 



P-l 0-1 1 (x 2 +y 2 +2xyp xy ) 

— —= jj £ A x p S(x - X p ) £ A y q S(y - Y q ) e 2 dxdy (9) 



jj ^-p-v- -py / , 

2TT II — Pxy -oo p = 0 q = 0 



Thereafter, it is easy to compute the integrals over the x and y domains by evaluating them at the 
points where Kronecker's delta does not vanish. Therefore, Equation (9) becomes Equation (10): 

aD . P-lQ-l l(Xj3 + Yj* +2X p Yq p xy ) 

6R xy _ i VV A , A y e -2 ( i- Pxy) 



dp. 



xy 



2lt ll-pJyP-Ofl.O 



ZZw 1 do) 



and finally, the integration over the p xy domain is performed from -1 to 1. 

P-lQ- 1 iCX^+Y^+ZXpY^y) 



R, 



= ( 1 V V a* A y e 2 V-Pw) do 

xy J I Z-i Z-i p q y (11) 



-l 2tt Jl -p2 yP =o q=o 



From Equation (11) it is straightforward to include the impulse response of the system by 
considering several values of p xy for several delays (t): 
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Pxyij) p _ x Q-x l(x^+Y^ + 2X p YgP yy (T)) 

R xy (r) = j 1 ^^A*A^e 2 MyW) dp xy (T) (12 ) 

0 2TT II - Pxy(j) P = 0 q=0 

In some particular and simplified cases, Equation (12) can be obtained analytically, for instance in 
the case of quantifying with two levels (one bit, Q = P = 1 and Aq = Aq = 2 ). In this case, 
Equation (12) becomes the well-known solution stated in Equation (13) [9]: 

R xy (PxyO)) = ^arcsin (p xy (r)) (13) 

Otherwise, when the quantization scheme is more complex Equation (12) can only be analyzed 
numerically, such as having different numbers of bits, when signals are clipped, when the quantization 
steps are irregular, when the ADC has a non-linear circuit before that exhibits gain compression, when 
there is no symmetry between the positive and negative parts or many other possibilities. 

Figure 2 has been obtained from Equation (12) for several quantization schemes. Two q[-] functions 
are plotted for reference, the coarsest quantization scheme (i.e., 2 levels/1 bit, Equation (13)), and the 
ideal one (i.e., infinite levels, R xy (r) = p xy (r)). Three more functions are shown considering an ADC 
span window of 5a x>y (Vadc = 5cr A;V ), 15 quantization levels or 3.9 bits (log 2 15 = 3.9, note that it is not 
necessary to have an integer number of bits), but changing the non-linear function. Several g(x) 
functions have been considered such as equally spaced levels, using gain compression modeled by an 
hyperbolic tangent (g(x) = tanh (x)) or having randomly spaced levels. All results are different from 
the others, showing the need of a deeper analysis. 

Figure 2. Relationship between the non-linear and the ideal correlation for different 
digitization schemes. 
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The root mean square error (RMSE) coefficient between R xy (r) and p xy (j) has been chosen to 
measure the degree of linearity of q[-]. So, the RMSE has been selected to be the figure of merit for 
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the comparison of several quantization schemes. Obviously, the linearity of any q[-] function is better 
as the RMSE is closer to zero. For a given application, there is a maximum distortion of the RMSE that 
can be afforded by the q[-], and that is called MaxRMSE- If the relationship between the quantized and 
the ideal correlations is linear enough (RMSE < MaxRMSE), then the value of the quantized correlation 
matches with the ideal correlation and there is no need of using the function. Otherwise, if it is 

not linear enough (RMSE > Max RMS E) then, the obtained correlation has to be modified by the 
q~* [•] function in order to retrieve the real correlation. 

Figure 3 presents the RMSE computed using Equation (12), with different quantization levels and 
Vadc- Each curve has a minimum that corresponds to the optimum configuration of Vadc with respect 
a x , y . As the number of quantization levels increases, the minimum value of RMSE curve gets closer 
to 0, and the Vadc/Oxj range where the curves remain close to 0 increases as well, since the function 
q[-] is linear over a wider range of input powers. By inspecting Figure 3 it is clear that for V A dc < 2&x, y , 
clipping has a dominant impact over the non-linear correlation. This effect has the maximum impact 
when the clipping reduces the whole ADC to a 2 level decision (1 bit). As the clipping effect increases, 
the RMSE value converges to the RMSE of the 1 bit/2 levels quantization (RMSE =14.2%), which is 
the most non-linear quantization scheme possible. 

Figure 3. Root mean square error and ADC span window relationship for different equally 
spaced quantification levels. 
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On the other hand, plots also have a common trend in the region above the minimum, where the 
ratio Vadc/<7a,} increases. As the number of quantization levels increases, the ideal and non-linear 
relationship (<?[-])becomes more linear over a wider range of input power (the region around the 
minimum widens). But when Vadc/o"*,? , the RMSE asymptotically increases again up 14.2% (the 
value of the 1 bit/2 levels) for all the quantization schemes. This effect can be understood as fewer bits 
are being used as the VADc/o^y ratio increases, and it can be understood as a decreasing of the effective 
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number of bits. The maximum effect occurs when the Vadc is so spread over er T y that only two levels 
are effectively used for quantization. 

Another conclusion that can be drawn from Figure 3, it is that each quantization level has its own 
optimum ADC span window (Vadc) with relation of the standard deviation of the input signal. This is 
a critical design parameter and it has been summarized on Table 1 for some quantization schemes. 

Table 1 shows that the two-level quantization scheme is not sensitive to the Yadc^xj relationship 
and the RSME is always 14.2% and, on the other hand for 31 levels the relationship between quantized 
and ideal correlation is almost linear. Furthermore, as the number of levels increases the minimum 
RMSE decreases exponentially (Figure 4) and it occurs at a higher \ Av>d<yx, y ratio. 

Table 1. Summary of the critical design parameter VADC/ax,y. Relationship between its 
optimal configuration and the minimum root mean square error obtained for some 
quantization schemes. 



Levels 


V A Dc/a x ,y optimum 


RMSE minimum | % ] 


2 (1 bit) 


no optimum 


14.2 


3 (1.6 bits) 


0.7 


4.7 


7 (2.8 bits) 


1.8 


1.4 


15 (3.90 bits) 


2.2 


0.3 


31 (4.95 bits) 


3.0 


0.02 


63 (5.97 bits) 


3.5 


0.0066 



Figure 4. Minimum root mean square errors for different quantization levels, it decreases 
following an exponential trend, (a) RMSE decreasing vs. the quantization levels, and 
(b) RMSE in a semi-log axis decreasing vs. the number of bits. 
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Figure 5 shows the effect of using various g(x) functions. It compares an equally spaced 
quantification with a hyperbolic tangent spacing, used to compress the gain. 
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Figure 5. Root mean square error for different levels and different spaced quantification 
levels, (a) 7 quantization levels, and (b) 15 quantization levels. 



RMSE for 7 quantization levels 

























































Equally 

Gain cc 


spacec 
m press 


ion 







































































































10 12 14 16 18 



20 



Wx,y 

(a) 



0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



RMSE for 15 quantization levels 















































































































































Equally spaced 

Gain compression 

















10 12 14 16 18 



20 



(b) 



In both cases, Figure 5(a) (seven levels) and Figure 5(b) (15 levels), there are no changes in the 
clipping area, but there are as the Vadc/Gxj ratio increases. The effect is that an ADC with gain 
compression has lower RMSE values than the same ADC but with equally spaced level scheme for a 
wider area. This is equivalent to increasing the number of quantization levels. An ADC with gain 
compression can be useful to ensure linearity over a wide range of input signal power, reducing the 
necessary number of bits. A common radiometric application where a high range of input power is 
measured is during the radiometric T hot ~ T co i d calibration. 

4. Impact of Quantization on the Correlation Spectrum 

In [7] some analytically easy-to-solve quantization schemes were analyzed using Taylor series 
approximations when p xy «l, which applies in most radio-astronomical measurements, but not 
necessarily in Earth remote sensing or during a radiometric calibration. The cross-correlation 
spectrum, the frequency power density is spread and distorted due to quantization step and clipping. In 
this work, the spectrum is computed following Equation (12). Therefore, the analysis can be 
exhaustive and include different quantization levels for two input signals for any complexity, gain 
compression, different bandwidths, and different channel frequency responses. The frequency power 
density is the Fourier transformation of the cross-correlation function, which is the Wiener-Khinchin 
theorem (Equation (14)): 

oo 

Sxy(f)= j R X y{T)e-^fdT = F{ gi {x{t))}F{g 2 {x{t))Y (14) 

— oo 

Furthermore, Equation (14) also shows the relationship of the cross-correlation spectrum with the 
individual signals spectrum, which is the multiplication of the Fourier transforms of the output signals 
of the quantization functions. By replacing Equation (4) in Equation (14) the spectral power density of 
the non-linear cross-correlation can be readily obtained (Equation (15)): 
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S xy (f)= j q[p xy {T)]e-^f dx (15) 

— oo 

where the function q[p xy (j)] can be obtained analytically (following Equation (13)) or numerically 
(following Equation (12)). 

The effect of the quantization in the cross-correlation spectrum is to decrease and distort it in the 
pass-band and spread it over the rejected band. Figure 6 presents an example of the impact of the 
quantization on the spectrum computed following Equation (15). The input signals, x and y, that have 
been considered for this analysis have a rectangular frequency channel's response with a pass-band 
bandwidth of 2B and a Vadc = 5er Y , v relationship for all the quantization schemes. 

Figure 6(a) presents the x and y linear cross-correlation on the delay domain for several quantization 
levels. The lower the number of levels, the higher the correlation distortion. Figure 6(b) shows the 
corresponding spectrum of cross-correlations presented in Figure 6(a), which have been computed 
following Equation (15). As expected, the two-level scheme has the highest spectrum distortion, the 
spectrum spread increases on the rejected band and the pass-band decreases. 

Figure 6. Effect of the quantization on the cross-correlation spectrum, (a) ideal 
cross -correlation (p xy (j)) and cross-correlation for different quantization levels, and 
(b) their corresponding spectrum. 



FWF for V ADC =5a x FFT(Cross-Correlation) for V adc =5ct x 




Delay Frequency 

(a) (b) 

Table 2. Summary of the spectrum distortion and spread for V A dc = 5c XJ and different 
quantization schemes. 



Levels 


Pass-band 


Rejected band increase 


2(1 bit) 


80% 


13% 


7 (2.8 bits) 


87% 


5.0% 


15 (3.90 bits) 


97% 


1.3% 


31 (4.95 bits) 


99% 


1.0% 


Ideal 


100% 


0.0% 
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Table 2 summarizes these results, the main conclusion is that using Vadc = 5er VjV and at least 15 
quantization levels it can be assumed that it neither spreads nor distorts the spectrum. Figure 7 shows 
the effect for two g{x) functions and for seven quantization levels. It can be observed that the 
compression gain reduces the distortion and the spread spectrum. 

Figure 7. Effect of the quantization on the cross-correlation spectrum, (a) ideal 
cross -correlation ( Pry CO) and cross-correlation for two different quantization schemes 
using 7 levels, and (b) presents their corresponding spectra. 



FWF for V ADC =5a x for 7 lewis FFT(Cross-Correlation) for V ADC =5o x for 7 levels 




-3/(2B) -1/B -1/(2B) 0 1/(2B) 1/B 3/(2B) -4B -3B -2B -B 0 B 2B 3B 4B 



Delay Frequency 

(a) (b) 



In this section, it has been shown that the quantization can be successfully expressed and analyzed 
as a non-linear transformation. Despite the quantization can affect significantly the value of the 
cross-correlation, it is possible to recover the ideal correlation using the inverse function which 
can be obtained numerically. Furthermore, for each quantization scheme there is an optimum 
configuration, in terms of linearity, of the Yadc/^x^ relationship. Furthermore, quantization distorts and 
spreads the cross-correlation spectrum. Finally, the quantization not only has an impact on the value of 
the cross-correlation and on the spread of the spectrum, but it increases the standard deviation of the 
cross-correlation, as it will be shown in the next section. 

5. Sampling Impact 

The main effect of sampling is that it creates replicas of the spectrum. Usually, in a standard digital 
radiometer topology there is an anti-aliasing filter before the digitization process. Even in this case, the 
digitization can introduce aliasing due to the spread of the quantized spectrum, which depends on the 
sampling frequency (F s ), and the baseband bandwidth (B). To understand how this effect impacts on 
the sampling it is necessary to define the Signal to Noise Ratio (SNR) as in Equation (16): 

SNR = -^ (16) 
where R xy is the expected non-linear cross-correlation value, and 0/? is the standard deviation. 
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Moreover, the SNR term is related with the radiometric resolution [18], which is one of the figures 
of merit of any radiometer. The SNR can decrease due to an increase of the standard deviation of the 
correlation or due to a decrease of the Rxy> i.e., a noise increment on the pass-band (Equation (2)). 

Figure 8 shows an ideal baseband spectrum with an infinite sampling frequency (F{p X y(j)}, in 
black) and a cross-correlation spectrum spread due to quantization (F[R xy (r)}, in green) for different 
F s and its replicas [19,20]. After quantization and sampling, and further cross-correlation of the 
signals, the quantization tail of the Nyquist's replica (F s = 2B, in red) overlaps the pass-band of the 
baseband replica, so an extra amount of noise is included in the measurement of R xy , and the effect is 
a decrease of the calculated mean value (Equation (2)). On the other hand, the standard deviation remains 
constant, but the SNR decreases, which has a large impact when considering sub-Nyquist sampling. 

Figure 8. Effect of the sampling on the correlation spectrum, (F{-} stands for Fourier 
transform). Spectrum of the ideal and the quantized correlation at baseband, black and 
green, respectively. In red, it is the first spectrum replica for the minimum sampling 
frequency which fulfills the Nyquist criterion (Fs = 2B). In blue, it is the first spectrum 
replica for Fs » 2B. 



F{Rx y (t)} 



Xoise increment 



F{px T (t)} 
/ Fs-2B replica 

■A / 



~Fs»2B replica 



Frequency 



If the sampling frequency increases, then the first replicas (F s » 2B, in blue) appear further away 
from the band-pass replica and the tail effect has less impact than before. This is why when only few 
quantization levels are used an increase of the sampling frequency has a significant impact on the 
radiometric resolution. Even though, it has less impact when more quantization levels are used. 

Figure 9 shows the impact on the quantized correlation of different sampling frequencies. The 
considerations for these plots are a 2 quantization levels scheme, a noise equivalent band-pass 
bandwidth of 2B (assuming a band-pass signal), and Vadc = 5er v , v , taking into account the previous 
considerations about the spectrum replicas. Figure 9(a) shows the deformation of the fringe washing 
function for a set of sampling rates. The sub-Nyquist sampling (F s = 0.75 FN yqu ist) has the largest 
deformation: it is wider, sharper than the rest of them, and the side lobes are higher. 

For the exactly Nyquist sampling rate the fringe-washing function is still sharp and wide. For a 
higher sampling rate (F s > 2 FN yqu ist) the function is not affected by the sampling rate. On the other 
hand, the spectrum in Figure 9(b) is totally distorted in the sub Nyquist sampling rate case, and 
recovers its original shape as the sampling rate increases. 
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Figure 9. Impact on the quantized correlation of different sampling frequencies. 
For 2 quantization levels, a bandwidth 2B (assuming a band-pass signal), and Vadc = 5a x , y . 
(a) shows the fringe-washing function deformation, and (b) shows the effect on their spectrum. 
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Figure 10 shows the impact on the quantized correlation of different sampling frequencies. The 
difference from the previous analysis is that in these plots 31 quantization levels have been considered. 
As this quantization scheme spread less the spectrum than the previous one, a better performance is 
achieved for lower sampling rates. Figure 10(a) shows the distortion of the fringe-washing function 
considering different sampling rates. The sub-Nyquist sampling (F s = 0.75 F Nyqu i st ) has the largest 
distortion: it is wider, sharper than the rest of them, and the side lobes are higher. For a larger sampling 
rate (F s > FN yqu ist) the function is not affected by the sampling rate. On the other hand, in the 
Figure 10(b), the spectrum is completely distorted for the sub-Nyquist sampling rate case, and recovers 
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its original shape as the sampling rate increases. As expected, it has better performance for higher 
sampling rates. 



Figure 10. Impact on the quantized correlation of different sampling frequencies. 
For 31 quantization levels, a bandwidth 2B, and Vadc = 5a XJ . (a) shows the deformation 
for the fringe-washing function, and (b) shows the effect on their spectrum. 
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6. Sampling Rate Inaccuracies 

In the conversion from analog to digital signals, the sampling frequency is normally assumed to be 
constant. Nevertheless, the sampling clock is prone to interferences, thermal drifts and other effects 
that can change the oscillating frequency and introduce sampling inaccuracies. 
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Figure 1 1 shows the effect of these inaccuracies, that can be summarized into skew (Tskew), which is 
the delay offset between the sampling time of x and y, and jitter, which is the time fluctuation above 
the mean value (T jitter). The impact of the skew is changing the point where the FWF is evaluated and a 
phase rotation of the correlation coefficient if the working frequency is not exactly zero. Meanwhile, 
jitter creates a perturbation for each correlation sample around the mean point with zero mean. 

Figure 11. Impact of the clock inaccuracies on the correlation value. T s k ew is the sampling 
rate offset between x and y, and Tji tt er is its fluctuations due to clock inaccuracies. 
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A clock skew error not only has an effect on the modulus of the correlation coefficient, but it affects 
its phase by rotating the real and imaginary parts [21], as it is shown in Equation (17): 

<x(nT 1 )y*(nT 2 )> 

= r xy {\ - T 2 ){[p R coso) r (T 1 - T 2 ) - p I smui r (T 1 - T 2 )] ^ 



+ y'[p / coscri r (T 1 - T 2 ) + pRsinUrCTi - T 2 )]} 



where: 

• (xinT^y* (nT 2 )> = p = p R + jp,, the correlation coefficient is a complex value with modulus 
and phase (j = V— T) 

• o) r = 27T-7 is the ratio between the residual frequency (f r ) and the sampling frequency, 
residual digital frequency of the signals, and F s = 1/T S 

• r xy(Ji ~ T 2 ) is the fringe-washing function evaluated at T x — T 2 , and 

• T x and T 2 are the sampling periods for the x and y inputs, respectively. Both sampling times 
are identical except for the skew and jitter effects ((T x - T 2 ) = T skew + Tjitter(O)- 
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On the other hand, the jitter is a random variable effect that can be statistically modeled as a 
polynomial of standard deviation depending on time (^jitter = °app + V c cik( n ~ m )Ts), where n and 
m are two different samples. So that, it can be associated to two different and independent 
processes ([22,23]): 

• the aperture jitter, which is related to the random sampling time variations in the ADC 
caused by thermal noise in the sample&hold circuit. The aperture jitter is commonly 
modeled as an independent Gaussian variable with zero mean and a standard deviation of 

^Jitter = tfapp* an d 

• the clock jitter, which is a parameter of the clock generator that drives the ADC with the 
clock signal. The clock jitter is modeled as a Wiener process, i.e., a continuous-time 
non- stationary random process with independent Gaussian increments with a standard 
deviation equal to Oj itter = ^c clk (n — m)T s . 

Taking into account a wideband signal in jitter terms (fulfill the inequality 27rBaj itter > 1). From 
here, x can be defined using its discrete-time Fourier transformation (Equation (18)): 

x(nT s + T /n ) = T s j Xtne-W^^JnMf (18) 

-F s /2 

where Ty n is the value of the instantaneous jitter. Thence, following Equation (18) the correlation can 
be defined again using Equations (14) and (18): 

F s /2 

R xy (r) = <x(mT s )y*(nT s + T /|m _ n| )> = ( j (x(f)Y*(f)e i2nfTj \™-n\) e -j2u/(m-n)T sd/) (19) 

-fs/2 

fs/2 

R xv {t)= j S xy {f){e ]2nn ^-n\)e-^ m - n ^Af (20) 

-F s /2 

where: 

• (e )2TT ^ T; i m - n i ) is the expected value of the jitter and it depends on the jitter type [23], 

• T/| m -n|i s the existing jitter between the m th sample of x and the n ,h sample y, and 

• t is the digital delay between the cross-correlation samples (t = (m — n)T s ). 

Equation (21) applies for an aperture jitter and Equation (22) applies for the clock jitter: 

( e J' 2TI / T 7a PPn ) = e 2n 2 f 2 a^ pp> md ^i 2 ^ JapP\ m - n} ) = g-2Ti 2 / 2 2^ pP) j f ^ dse j (2 1) 
( e j2Tt/T ;clkn ^ _ e 2n 2 / 2 c c i k nT S)and ^ e ) 2TT -^ T ycik| m _ n ^ _ e -2n 2 f 2 c clk T s \m-n\ Q2) 

where c cUi is the phase noise constant of the oscillator. Thence, the impact of the aperture jitter to the 
cross -correlation coefficient is shown in Equation (24): 



R xy (r)= j 5 %y (/)e- 2 - 2 ^ 2 ^te re -j2u/ M/ 

— oo 

oo 

= j S xy (f)e~ 2n2f2 ^ 2 ~ Sx)a ^pp +CclkTslm ~ n ^ e~^ 2nfz df 



(23) 



— oo 
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e 2a a PP e 2c clk n 
R xv (t) = R xv (t) * * - ■ 



^]2na app ^j2nc clk T 



where: 



• R xy (t) is the cross-correlation coefficient modified by the jitter, 

• 11 xy (T) is the cross-correlation coefficient, and 

• * denotes the convolution operator, 

From Equation (22) it can be inferred that the aperture jitter does not depend on the sampling rate 
and it has a low-pass filter effect on the correlation spectrum with its cut-off frequency depending on 
^app- Particularizing Equation (24) in the case of having a rectangular spectrum, then Equation (24) 
becomes Equation (25), and the cross-correlation function is multiplied by an erf(x) function 
(erf(x)=^j; e - t2 dt): 

B/2 



R xy (r) = a x a y p(T)e-^ 2 ^ 2 - s ^ppe-^df = R xy (t) 

_{ /2 4 2na *vv 



(25) 



As it can be seen, in the case of aperture jitter the coherence loss ratio (Equation (20)) is 
independent of the number of sampling points N (therefore of the integration block length TyV), while 
in the case of clock jitter it strongly depends on it (Equation (21)). 

7. Analysis of the Cross- Correlation Variance due to the Digitization 

In the previous sections, quantization and sampling have been analyzed in terms of the impact on 
the retrieved value of the correlation. In this section, the variances are analyzed taking into account 
both effects (quantization and sampling). In a general case, the output of a correlator of two quantized 
and sampled signals after averaging N q samples is defined as: 

N q 

R X y(t) =yY j 91 ^ x( - nT ^9Ky(mT s )) (26) 



where: 



• T s is the sampling period for the x and y inputs, respectively. In most of the 
applications both sampling times are identical except for skew and jitter effects 
«Ti - T 2 ) — Tskew + Tji tt er(0), for simplicity, without loss of generality, skew and jitter 
can be simulated as well), and 

• 9i,2 are defined in Equation (3). 

Obviously, the product of a pair of samples of P\ m - n \ = gi(x(nT s ))g2(y(mT s ) in Equation (26) 
does not follow a Gaussian pdf. Rather, the distribution of V\m-n\ follows a Rayleigh pdf when the g 
functions are linear (no quantization). In a more general case, having non-linear g 12 functions, the 
distribution of V\m-n\ follows different trends depending on the quantization scheme. However, the 
average of many of these products, over a large number of pairs of samples, approaches a Gaussian 
pdf, as implied by the central limit theorem. The variance of the correlator output can then be 
calculated as in Equation (27): 
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a R xy ~ (RxyRxy) (Rxy) 2 



(27) 



where the first term of this addition is expanded as Equation (28): 

N q N q 

(R xy R* xy ) = ^ ^ (9i{Knl s ))gl{y{ml s ))gl{x{nl s ))g 2 {y{ml s ))) (28) 

1 n =l m=l 

Thence, Equation (28) is split in two parts, the first one takes into account only the elements which 
have the same n and m time indexes, and a second part, which is the rest of the Equation (28): 

N q 

2<5i(*(nT s ))flf5(y(nT s ))flf;(x(nT s ))flf 2 (y(nT s ))> + 



(R xy R xy ) — 



n=l 



1 1 

^ ^ <5i(*(nT s ))flf5(y(nT s ))flf;(x(mT s ))flf 2 (y(mT s ))> 



n=l m=l 



(29) 



the first term of Equation (29) can be computed using the properties of the fourth and second statistical 
moments relationship for zero mean Gaussian random variables [18], as follows: 



{g 1 (x(nT s ))g* 2 (y(nT s ))gl(x(nT s ))g 2 (y(nT s ))) = 

= f <gdx(nT s ))gXx(nT s )))(gZ(yCnT s ))g 2 (y(7CT s ))) ] = 
[+(9i {x{nl s ))g* 2 {y{nl s ))){g{ {x{nl s ))g 2 (y (nT s ))> J 

= °L(x)°g 2 (y) + UiPxy(0)]\ 2 ^ l(x) ^ 2(y) = Og^CTg^ (l + |q[Pxy(0)]| 2 ) 



(30) 



where: 



• o^ x (x) an ^ <7 g 2 {y) are tne var i ances of 9i 0*0 an( l 9z(y)> respectively, and 

is the modulus of the cross-correlation coefficient of the input signals after the 
non-linear function. 

The second term in Equation (30) can be expanded using the same statistical moments relationship: 

N q N q 

a^Z Z^ l(x(nTs)) ^ (y(nTs)) ^ (x(mTs)) ^ 2(y(mTs))> = 



1 n=l m=l 



N n -1 



^y{g 1 (x(n%))gKx(m%)))(gKy(n%))g 2 (y(mT: s ))) 

N « rf=l 

Nq Nq 

+ A^Z Z ^^ x( ; nT ^)^(y(nT s ))>(^(x(mT s ))^ 2 (y(mT s ))> 

^ n=l m=l 

= 1 ^ L ~ kk(°)]| 2 <W<(y) + ^Z^ " n) ^( nT s) R yy(" T s) 



(31) 



where: 



R xx (riT s ) and i? yy (nT s ) are the auto-correlation of the g^x) and g 2 {y), respectively taking 
into account all the samples between sample 1 and the N q , as in Equation (12). 
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The first part of the variance of the correlation is obtained by combining the results obtained in 
Equations (30) and (31): 

-(l + |q[p*y(0)]| ) + 

N q 

, r „ — l 2 2 

1 9 [Pxy (°)] I o£ w o| 2 (y) + 2j ^ ~ n ^ Rx% (nT ^ R yy (nT ^ 



°gi(*)°g 2(y) 



n=l 



°"gi W°gz(y) 



2 x -1 2 

+ ^2 2j ~ n ) ( nT s) R yy ( nT s) + 1 9 [Pxy (0)] | o£ (x) C7^ 



2 

'52 (y)- 



n=l 



(32) 



Finally, the value of the variance of the N q -samples averaged correlation is given by: 

a 2 a 2 

<y = + ^2 X ^ - ^ RxxWRyyW 



(33) 



and arranging Equation (33), the variance appears in its usual form, as it is known in the literature [17]: 



a ly = F" + W T ( N ct - n) R xx (.m s )R yy (n%) 



(34) 



where, — Og x (x) a ]j (y) ls me var i ance f° r eac h sample of the correlation (N q = 1). 

2 

Following Equation (25), Figure 12 shows the evolution of °R xy with the integration time. The 
signals considered in these plots have a noise equivalent squared bandwidth of pass-band of IB, 
Pxy(°) = 1> and v adc = 5cr x , y . 

Figure 12. Impact of the number of averaged samples on the ctr , (a) shows the impact 
for a 3 level quantization scheme and different sampling rate, and (b) shows a constant 
sampling rate changing the number of quantization levels. 
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Figure 12. Cont. 
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Figure 12(a) shows the results for the three quantization levels, all the curves converge to zero 
asymptotically. The F s = 0.75FN yqu ist curve converges faster, because its samples are uncorrelated 
among them. For a given number of averaged samples, the variance of the quantized correlation 
increases with the sampling rate, because successive samples are more and more correlated. Recall 
that, abscises axis of these plots are displayed in number of averaged samples, not in integration time 
units. Obviously, for a gi ven integration time, the plot with the highest sample rate has the 
lowest a R x it has had time to average more samples despite that for a given number of samples it 
exhibits the worst performance. 

In Figure 12(b) the sampling rate has been frozen to F s = F Nyqu i st and two quantization schemes has 
been analyzed (3 and 15 levels). The three-level curve converges more quickly to zero than 
the 15 levels curve. This is because the fringe washing function is shaper due to the higher non-linear 
distortion in three levels than in the 15 levels quantization. This effect is equivalent to have more 
uncorrelated samples. Despite the mean of the correlation is lower for the three-level that the 15-level 
quantization, the variance of the correlation is also lower in this case. 

8. Conclusions 

The impact of a general quantization scheme has been expressed and analyzed as a non-linear 
transformation. Despite the quantization can significantly affect the value of the ideal cross-correlation 
(Pxy), it is possible to recover the ideal correlation (p xy ) from the measured one (R xy ) using then the 
function p X y = <? _1 [^%y], which can be always numerically obtained. It has been found that for each 
quantization scheme there is an optimum configuration of the Wadc/pxj relationship for which the 
linearity with respect the ideal correlation is minimized. Furthermore, quantization distorts and spreads 
the cross-correlation spectrum, decreasing the radiometric resolution. 
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Sampling also has an impact in the spectrum of the cross-correlation function, and depending on the 
sampling rate, replicas of the spectrum can overlap the main spectrum. This can be due to a sampling 
rate below the Nyquist criterion or, even above the Nyquist criterion, due to the non-linear quantization 
process which spreads the spectrum. In both cases, there is an impact on the SNR, and so on the 
radiometric resolution. Other sampling effects have been analyzed such as the clock inaccuracies: 
skew and jitter. 

Finally, the impact of quantization and sampling on the variance of the measurements has been 
studied. It has been found that it strongly depends on the relationship between the bandwidth of the 
cross-correlation and the sampling rate. As the sampling rate increases, the successive samples are 
more correlated, so they add less new information to the measurements and the variance 
decreases slowly. 
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